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We consider the limit of strong Coulomb attraction for generalized 
Hubbard models with ^-symmetry. In this limit these models are 
equivalent to the ferromagnetic spin-1/2 Heisenberg quantum spin 
chain. In order to study the behaviour of the superconducting phase in 
the electronic model under perturbations which break the ^-symmetry 
we investigate the ground state of the ferromagnetic non-critical XXZ- 
chain in the sector with fixed magnetization. It turns out to be a large 
bound state of N magnons. We find that the perturbations considered 
here lead to the disappearance of the off-diagonal longe-range order. 
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There has been a lot of activity in the past few years in the investigation 
of electronic models for superconductivity. Especially the so-called ^-pairing 
mechanism introduced by Yang ]I| has been quite fruitful for the construction 
of models with superconducting ground states (see e.g. @). Although the 
ground state of the Hubbard model is not of the 77-pairing type |J one can 
construct such models by adding nearest-neighbour interaction terms. The 
first example was found in || |J . In [Q] it has been shown that for a large class 
of generalized Hubbard models the ground state in the limit U — > —00 is the 
simplest 77-pairing state {rf) N \fy with rf = J2j c ]j. c jT' F° r the supersymmetric 
Hubbard model |3], || and a special case of the Hirsch model with correlated 
hopping interaction |j) the complete phase diagram in arbitrary dimensions 
has been found. It shows two phases with 77-pairing ground states. Since 
these states have ODLRO they are also superconducting ||. 

The main aim of this letter is to clarify the effect of //-symmetry-breaking 
perturbations on the superconducting ground state. In order to keep the 
problem as simple as possible we restrict ourselves to the one- dimensional 
case and the limit U — > — 00. For strong coupling limits spin models often 
turn out to be effective Hamiltonians for models of correlated electrons. E.g., 
the antiferromagnetic Heisenberg model is known to be the effective Hamil- 
tonian for the Hubbard model in the limit U — >• 00 at half-filling. In our case 
the effective Hamiltonian is found to be the ferromagnetic Heisenberg chain. 

We perturbate the supersymmetric Hubbard model by changing the value 
of the nearest-neighbour Coulomb interaction V in such a way that the ef- 
fective Hamiltonian in the limit of large attraction will be a ferromagnetic 
non-critical XX Z model (A > 1, see eq. (§) below). In this regime the 
ground state is simply the fully polarized ferromagnetic state \\ ■ ■ ■ f ), an d 
the spectrum has a gap which vanishes in the limit A — > 1 + . The ques- 
tion of stability of the superconducting ground state of the electronic model 
thus turns out to be closely related to the form of the ground state of the 
ferromagnetic XXZ model at fixed magnetization. 

In the following we will be interested in the large-?/ limit of generalized 
Hubbard models with Hamiltonian 
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Here c] CT and Cj„ are the usual electron creation and annihilation operators, 
i.e. {cj a ,cl a ,} = 5ji5 a(7 i, rij a = Cj a Cj a is the corresponding number operator 
and the components of the spin operators are given by S* = | (n^ — n,j^), 

S? = h (St + St) and S? = I ( St - S~) with St = clc,, , Sr J 



. = - ^ . + b . j ana t>; = s (^ - ^ j witn ^ = c^, ^ = C ^c iT . 
In [|J it has been investigated under which conditions the Hamiltonian 
([!]) has an ^-pairing ground state of the form 
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withr^t = Ef =1 v] = ZU4Av 

t = X,2V = Y,V < and 


It has been found that this is true for 



-[/>2max(V-^, V+^ + \lf\, 2V + 2|t|). (3) 

Note that the first two conditions also guarantee the //-symmetry of the 
Hamiltonian (D with U = 0, i.e. [H(U = Q),rj*] = 0. 

For the supersymmetric Hubbard model of f||, |4| only the Coulomb in- 
teraction U is a free parameter. The other interaction constants have fixed 
values: X = t, V = —t/2, Y = —t, and J xy = J z = It. The other inte- 
grable model ([!]) has - apart from the Coulomb interaction U - only two 
non- vanishing interactions constants, namely t = X ||. 

In the limit of strong attraction (U — > — oo) in the Hamiltonian ([!]) one 
expects [] a gap for the single-particle excitations which is proportional to \U\. 

^his can be shown explicitly for the exactly solvable cases M, H. 



This means that only the motion and interactions of the doubly occupied sites 
are important. Therefore all the terms in (JJ) can be dropped except for the 
pair-hopping term Y and nearest-neighbour Coulomb interaction V. Setting 
Y = — 1 (since the ^-symmetry implies Y < 0) we then get the effective 
Hamiltonian 
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(4) 
with A = 2V/Y = 1. Allowing for an interaction constant A/ 1 thus means 
a perturbation of the original model. In this letter we want to investigate 
the effects of A > 1 (i.e. 2V > Y). Note that this perturbation destroys the 
^-symmetry of the original Hamiltonian. 

In order to derive an effective spin Hamiltonian we first make a partial 
particle-hole transformation c-^ — > cL and cj, — > Cj± on the j-spins only. 
This changes n^ — > n^ and riji — > 1 — riji and transforms doubly-occupied 
sites into j-spins and empty sites into j-spins. The transformed effective 
Hamiltonian (f|) now reads 

+ c jT c il- c i+U c ?+ 1 >T + c ii c -J't c i+i,T c i+ 1 >lJ- W 

Using the spin operators introduced above and S 1 " = |<r" (where cr" are the 
Pauli matrices, a = x,y, z) the Hamiltonian (|5|) becomes 

^e// = -^ E (<W+1 + Ojof+l + M°i+l) ( 6 ) 



Under the particle-hole transformation - which maps the fermionic vac- 
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uum | ) onto the ferromagnetic state | J, ■ • • j) = Il,=i cl, | ) - the ?7-pairing 



state (0) translates into 



M = (s + ) N \i-i) (7) 



which is the ground state of the isotropic Heisenberg ferromagnet in the 
sector with fixed magnetization S z = N — ~. 



In the following we will determine how the form ([?]) of the ground state 
changes for A > 1 by proving a conjecture by Gaudin ||. 

Hamiltonian (||) has long been known to be solvable by means of Bethe 



Ansatz || |10J. In the sector with a fixed number N of upturned spins, the 
eigenvalues are determined by the equations [[U], [11], |12| 
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We will consider all possible types of string solutions 

AK = A^-^(nH-l-2fc) (A; = 1, 2, . . . ,n with n E Z+). (10) 

Periodicity along the real axis allows one to take the real center of the string, 
A£t\ in the restricted range (— n/2, tt/2]. As usual, it is possible to reduce 
(H) to a set of equations for the string centers. Choosing 
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we find 
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where M m denotes the number of m-strings, /^ is integer (half-odd) if L + 
M n + 1 is even (odd) and 
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The energy of a n-string, derived from (P), is 

e „(A)=2sinh 7 **"" - (14) 

cosnrry — cos 2 A 

and therefore the total energy, neglecting the immaterial additive constant, 
is 

E = 2sinh 7 V V ^^ r- (15) 

n=1 p =1 cosh Try — cos 2A^ 

where {A« } are the real string centers. We are interested in the ground 
state in a sector of fixed magnetization where 
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In principle one should find all solutions of the Bethe Ansatz equations for the 
string centers A^ n ) which are compatible with (|1^) and compare the respective 
energies. Yet, for a string of length n, lower and upper bounds of the energy 
( |HD are given by 

. sinhrry . sinhrry . sinhrry 

2 sum 7 < 2smli7 < 2smh7 (17) 

cosh wy + 1 cosh 727 — cos 2A cosh 727 — 1 

for A = 7r/2 and A = respectively. We conclude that the ground state must 
contain exactly one string of length iV if we can show that, for any other 
string configuration (Mi, M 2 , . . . Mjv-i, 0) with Yln=i n M n = N, we have 

sinh A^7 1 lz^ , , sinh wy . , 

A7 - < Y M n '— (18) 

cosh Nj — 1 ~^ cosh wy + 1 

First we observe that, in the sector under consideration, any configuration 
(Mi, M2, . . . , Mjv-i, 0) can be obtained starting from one of the configura- 
tions where two strings only are present, M n = 1, M m — 1, n + m — N, and 
performing steps in which a string is broken into two smaller strings 
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(Mi,...,M ni ,...,M n2 ,...,M n 

-»■ (Mi, . . . , M m + 1, . . . , M n2 + 1, . . . , M n - 1, . . .) (19) 



with rii + ri2 = n. In each step the RHS of ([H]) increases since 

sinhra 1 7 sinhn 2 7 sinhn7 . . . N 

H r r> ; 7 (ni + n 2 = n). (20) 



coshni7 + l coshn27 + l coshrry + l 

This inequality is guaranteed by the fact that /(0) = and f"(x) < in 
(0, +00), where f(x) = J^^+i > and consequently f(x 1 ) + f(x 2 ) > f(x 1 +x 2 ). 
We conclude that, to demonstrate flT8| ) it is sufficient to consider configura- 
tions on the RHS made up of two strings only. Again concavity of f(x) shows 
that it is enough to have 

sinhA^7 sinli7 sinh(iV — 1)7 



cosh A/7 — 1 cosh 7 + 1 cosh(iV — 1)7 + 1 

and, for fixed 7, this is certainly the case for TV large enoughP]. This concludes 
the proof that the ground state is given by a solution of the Bethe Ansatz 
equations with only one iV-string confirming a result predicted by Gaudin 
[§]. For the isotropic case A = 1 a similar result was already obtained by 
Bethe @. 

The string center A^ of the Bethe state with a single iV-string flT0| ) can 
easily be determined from fll2"l). On a finite lattice there are L solutions 

,(N) ( , ^7 kI (N) \ 
A v ; = arctan tanh tan , 
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(this holds for both L even and L odd) with the ground state at I^ N ' = -| and 
energy given by fll4|). The energies of all these states become degenerate in 
the limit L — > 00 (and consequently N — > 00) with limiting value 2 sinh7 and 
energy differences vanishing like 0(e~ Nl ). Consequently we have an infinite 
degeneracy of the ground state. 

Low-lying excitations are given by configurations where the iV-string is 
broken into a finite number of shorter strings. In this situation, the "in- 
teraction term" in equation (|12|) for the centers A^ goes to zero like 1/L. 
Therefore, the energy of each string coincides with the bare energy (|14D and 



2 Note that for large TV the difference "w^T^L - si ! lh J V7 -■, is of order e^l 

cosniiV — II7+I cosh A* 7 — 1 



no dressing needs to be considered. In the most general excited state con- 
figuration (Mi, M 2 , . . . , M/v-i) with Y<n=i nM n = N, the total number of 
strings Yln=i M n remains finite and the length of one or more strings has 
to diverge. Each of these diverging length strings contributes 2 sinh 7 to the 
energy, regardless of their position, while the finite length ones have energy 
fll^l). It is easily seen that the spectrum has a gap 

AE = 2(cosh7- 1) (23) 

obtained by taking Mi = 1, M N _ X = 1. 

The above results show that the ground state of the ferromagnetic XX Z 
chain (^j) with A > 1 at fixed magnetization S z = N — |f is a bound state of 
the form 

\M= £ #ci,...,^)n<.|l •••!)■ (24) 

{xx,...,x N } j=l 

Here ip(xi, . . . , xn) is a bound state wave function, i.e. it decays exponentially 
with respect to all coordinate differences \xj — xi\ — > 00 (J, I — 1, . . . , iV with 
j ^ I) |§. The wave function ip(xi, . . . , xjy) for x\ < x 2 < ■ ■ ■ < x^ is given 
by the Bethe-Ansatz expression 
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(25) 
where SV denotes the permutation group. We now specialize to the state 
containing exactly one iV-string. It can be seen easily that only the term 
where P is the identity has a non-vanishing contribution. Introducing the 
new variables z Q = -^ J2f=i x j an d Zj = Xj + i — Xj (j = 1, . . . , N — 1), i.e. 
Xj = zq + -^ J2k=i kzk — J2k=l z k we finally can rewrite the wave function as 
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where the constant C(A, zo, N) depends only on A, zq and N . The form ( p6|) 
shows directly that the wave function describes a bound state since it decays 
exponentially as function of the coordinate differences Zj = Xj + \ — Xj. 



For the fermionic model the ground state is then of the form 

N 

\1>n)= E V^i,...,a*)n(4,l4,T)l >- ( 27 ) 
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This state has no ODLRO, i.e. — J , l , : '\ ,' T ! l — — — ► 0. This can be 
seen as follows: 

Using the SU(2) commutation relations for r)j = Cj-jCj-j, || it is easy to see 
that 

where N is a normalization constant. Since ip(x\, . . . ,xn) decays expo- 
nentially as function of all coordinate differences a significant contribution 
to the sum in ( ffif) comes only from x\ ~ X2 ~ ... ~ Xjv-i ~ J and 
x\ ~ X2 ~ • • • ~ Xiv_i rj /. This means that for |j — Z| — >■ oo at least 
one of the factors in each term and thus the whole sum decays exponentially. 
This argument shows that superconductivity in the one- dimensional su- 
persymmetric Hubbard model is destroyed by a perturbation with A > 1. 
Instead, the form of the ground state suggests a tendency to phase separa- 
tion. This is not surprising since numerical investigations have shown that 
superconducting phases in electronic models appear quite generally in the 



vicinity of phase separation [13]. In higher dimensions however the situation 
is somewhat different. By analogy with the non-ideal Bose gas and prelim- 
inary results from perturbation theory [TJJ] we expect the superconducting 
ground state to be stable under the kind of perturbation considered here. 
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